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Abstract 

A new class of exact solutions of Einstein’s field equations with bulk 
viscous fluid for an LRS Bianchi type-I spacetime is obtained by using 
deceleration parameter as variable. The value of Hubble’s constant Ho is 
found to be less than one for non-flat model and is equal to 1.5 for flat 
model which are of the physical interest. Some physical and geometric 
properties of the models are also discussed. 
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1 Introduction 

The Bianchi cosmologies play an important role in theoretical cosmology and 
have been much studied since the 1960s. A Bianchi cosmology represents a spa¬ 
tially homogeneous universe, since by definition the spacetime admits a three- 
parameter group of isometries whose orbits are space-like hyper-surfaces. These 
models can be used to analyze aspects of the physical Universe which pertain 
to or which may be affected by anisotropy in the rate of expansion, for exam¬ 
ple , the cosmic microwave background radiation, nucleosynthesis in the early 
universe, and the question of the isotropization of the universe itself pQ. For 
simplification and description of the large scale behaviour of the actual universe, 
locally rotationally symmetric [henceforth referred as LRS] Bianchi I spacetime 
have widely studied m-m- When the Bianchi I spacetime expands equally in 
two spatial directions it is called locally rotationally symmetric. These kinds of 
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models are interesting because Lidsey [7] showed that they are equivalent to a 
flat (FRW) universe with a self-interesting scalar field and a free massless scalar 
Held, but produced no explicit example. Some explicit solutions were pointed 
out in references mm- 

The Einstein’s field equations are coupled system of high non-linear differ¬ 
ential equations and we seek physical solutions to the field equations for their 
applications in cosmology and astrophysics. In order to solve the field equations 
we normally assume a form for the matter content or that spacetime admits 
killing vector symmetries m- Solutions to the field equations may also be 
generated by applying a law of variation for Hubble’s parameter which was pro¬ 
posed by Berman m- In simple cases the Hubble law yields a constant value 
of deceleration parameter. It is worth observing that most of the well-known 
models of Einstein’s theory and Brans-Deke theory with curvature parameter 
k = 0, including inflationary models, are models with constant deceleration pa¬ 
rameter. In earlier literature cosmological models with a constant deceleration 
parameter have been studied by Berman HU. Berman and Gomide mi, Johri 
and Desikan mi, Singh and Desikan [14], Maharaj and Naidoo [15], Pradhan 
and et al. [TB] and others. But redshift magnitude test has had a chequered 
history. During the 1960s and the 1970s, it was used to draw very categorical 
conclusions. The deceleration parameter q 0 was then claimed to lie between 0 
and 1 and thus it was claimed that the universe is decelerating. Today’s situ¬ 
ation, we feel, is hardly different. Observations [mill] of Type la Supernovae 
(SNe) allow to probe the expansion history of the universe. The main conclusion 
of these observations is that the expansion of the universe is accelerating. So 
we can consider the cosmological models with variable deceleration parameter. 
The readers are advised to see the papers by Vishwakarma and Narlikar m 
and Virey et al. |20j and references therein for a review on the determination 
of the deceleration parameter from Supernovae data. 

Most cosmological models assume that the matter in the universe can be de¬ 
scribed by ‘dust’(a pressure-less distribution) or at best a perfect fluid. The role 
of the bulk viscosity in the cosmic evolution, especially as its early stages, seems 
to be significant. The general criterion for bulk viscosity was given by Israel 
and Vardalas [21], Klimek [22] and Weinberge ;,23j). For example, the existence 
of the bulk viscosity is equivalent to slow process of restoring equilibrium states 
(Landau and Lifshitz [24]). The presently observed high entropy per baryon 
ratio in the universe can be explained by involving some kind of dissipative 
mechanism (e. g., bulk viscosity). Bulk viscosity is associated with the GUT 
phase transition and string creation. Thus, we should consider the presence of 
a material distribution other than a perfect fluid to have realistic cosmological 
models (see Gr0n [25]) for a review on cosmological models with bulk viscosity). 
The model studied by Murphy [55] possessed an interesting feature in that the 
big bang type of singularity of infinite spacetime curvature does not occur to 
be a finite past. However, the relationship assumed by Murphy between the 
viscosity coefficient and the matter density is not acceptable at large density. 
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The effect of bulk viscosity on the cosmological evolution has been investigated 
by a number of authors in the framework of general theory of relativity (Pad- 
manabhan and Chitre m, Pavon [25], Johri and Sudarshan [25], Maartens [3D], 
Zimdahl m , Santos et al. 323, Pradhan, Sarayakar and Beesham [53], Kalyani 
and Singh [33], Singh, Beesham and Mbokazi [33], Pradhan et al. [3D]), Singh 
et al. 133 ), Bali and Pradhan [33]). This motivates to study cosmological bulk 
viscous fluid model. 

Recently Paul [35] have investigated LRS Bianchi type-I cosmological mod¬ 
els with a variable deceleration parameter. In this paper, we propose to find 
LRS Bianchi type-I cosmological models in presence of a bulk viscous fluid and 
we will generalize the solutions [39]. 


2 The Metric and Field Equations 

We consider the LRS Bianchi type-I metric in the form [5j 

ds 2 = dt 2 - A 2 dx 2 - B 2 (dy 2 + dz 2 ), (1) 

where A and B are functions of x and t. The stress energy-tensor in the presence 
of bulk stress has the form 


Tij ={p + p)uiUj - pg ij: 


( 2 ) 


where 

P = P~ (3) 

Here p, p , p and £ are the energy density, thermodynamical pressure, effective 
pressure and bulk viscous coefficient respectively and Ui is the four velocity 
vector satisfying the relations 

UiU 1 = 1, (4) 

The Einstein’s field equations (in gravitational units c = 1, G = 1) read as 

Rij —Rcjij 8nTij, (5) 


where R,j is the Ricci tensor; R = g l iR^ is the Ricci scalar. The Einstein’s 
held equations ([5]) for the line element <[T|) has been set up as 


2 B B 2 B' 2 

~B + B 2 ~ A2R2 “ ~ 8?rp ’ 


( 6 ) 


•, B'A 
B' -— = 0, 


A B AB 
A + B + AB 


A 

B" 


+ 


A'B' 


A 2 B A 3 B 


= -87 Tp, 


(7) 

( 8 ) 
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(9) 


2 B" 2 A'B' | B' 2 

~A 2 B ~ A 3 B + A 2 B 2 ’ 
The energy conservation equation yields 


P + (p + p) 



2 AB B 2 
AB ~ H 2 



—8irp. 


( 10 ) 


where dots and primes indicate partial differentiation with respect to t and x 
respectively. 


3 Solution of the Field Equations 

Equation 0, after integration, yields 


A = 


B' 

T’ 


(ii) 


where t is an arbitrary function of x. Equations © and © , with the use of Eq. 
m , reduces to 


B_d_ fB\ 


( 12 ) 


Since A and B are separable functions of x and t, so, is a function of x alone. 
Hence, after integrating Eq. m gives 

B = £S(t), (13) 


where S' is a scale factor which is an arbitrary function of t. Thus from Eqs. 
<ra> and m we have 

H = jS. (14) 

Now the metric 0 is reduced to the form 

ds 2 = dt 2 — S 2 [dX 2 + e 2X (dy 2 + dz 2 )] , (15) 


where X = In t. 
obtained as 


The mass-density, effective pressure and Ricci scalar are 


8n P = $2 [ 


8irp = 

R = 


1 

S 2 L 
_6 
S 2 L 


S 2 -l 
£2 


1 -S z - 2SS 
SS + S 2 - 1 


(16) 

(17) 

(18) 
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The function S(t) remains undetermined. To obtain its explicit dependence on 
t, one may have to introduce additional assumption. To achieve this, we assume 
the deceleration parameter to be variable, i.e. 


SS 


Q = ~ 


S 2 


(h + h 2 \ 

V H2 ) 


= b(variable), 


(19) 


where H = ^ is the Hubble parameter. The above equation may be rewritten 
as 

q q 2 

5 + t 5*-°- (20) 

4 Solution for b = — where a is constant 

In this case, on integrating, Eq. (l2()l) gives the exact solution 

S = iaf 2 + kt + d, (21) 

where k and d are constants of integration. Here we consider the two following 
cases: 


4.1 Non-flat models 

The mass-density, pressure and Ricci scalar are given by 


8n P = ^2 + at ) 2 “ 1 1 > 


8 t r(p -£8) = - 


/? = A 
R s 2 


i 

s 2 


( k T ut) 2 — 1 T 2n(— at 2 T kt -t- d) 


(k + at) 2 — 1 + a(-at 2 + kt + d) 


( 22 ) 

(23) 

(24) 


Here £, in general, is a function of time. The expression for kinematical param¬ 
eter expansion 6 is given by 


9 = 


3 (k + at) 

S 


(25) 


Thus, given £(t) we can solve the equations. In most of the investigations 
involving bulk viscosity is assumed to be a simple power function of the energy 
density 

m = Z0P n , (26) 

where £o a- n( l fi are constants. For small density, n may even be equal to unity 
as used in Murphy’s work for simplicity [7]. If n = 1, Eq. (1261) may correspond 
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to a radiative fluid [23]. However, more realistic models [32] are based on n 
lying in the regime 0 < n < i. 

On thermodynamical grounds, in conventional physics £ has to be positive; 
this being a consequence of the positive entropy change in irreversible processes. 
For simplicity and realistic models of physical importance, we consider the fol¬ 
lowing two subcases (n = 0,1): 


4.1.1 Model I: solution for £ = £o 

When n = 0, Eq. (l26l) reduces to £ = £o = constant. Hence in this case Eq. 
(l23l) . with the use of (l25l) . leads to 


8irp 


247Tg 0 

S 


(k + at) 


1 

S 2 


(k + at) 2 


1 + 2a(-at 2 + kt + d) 


(27) 


4.1.2 Model II: solution for £ = £o p 

When n = l, Eq. (l26l) reduces to £ = £o P- Hence in this case Eq. (l23l) . with the 
use of (USD, reduces to 


87 t p = 


{(k + at) 2 — 1} 

53 


[9^o (k + at) — S'] — 


2a(^at 2 + kt + d) 
S 2 ' 


(28) 


4.1.3 Physical behaviour of the models 

The effect of the bulk viscosity is to produce a change in perfect fluid and hence 
exhibit essential influence on the character of the solution. We also observe here 
that Murphy’s condition [7] about the absence of a big bang type singularity in 
the finite past in models with bulk viscous fluid, in general, not true. 

From Eqs. (l22l) and (l24l) . it is observed that 

1 — fc 

p > 0 and R > 0 for t > -, (29) 

where k < 1. From Eqs. (l22l) . (l24l) and (l25l) . it is also observed that p, R and 0 
decrease as t increases. 

We find that shear a = 0 in the models. Hence < 7 /9 = 0, which shows that 
the models are isotropic. From Hubble’s parameter equation, H = |[, we have 
an epoch time to given by 


to = 


_L * 

H 0 a 


which gives that 


at 0 + k = — -1 

Ho 


\J a 2 + H 0 2 (k 2 — 2ad) 
aH 0 


sja 2 + H 0 2 (k 2 - 2ad) 
Hn 


> 0. 


(30) 


(31) 
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From above equation we conclude that 


He 


- 1 > 0 , 


which reduces to 

H 0 < a. (32) 

From Eq. (13U1) . we observe that to > 0 for H 0 < | and k 2 > 2 ad i.e. 

Ho< Td' 

From Eq. ((Ml) . R > 0 implies that 

(k + at) 2 > [\J 1 — a(^at 2 + kt + d)) 2 . (33) 

It is evident from Eqs. (1M1) and (1331) that a < 1 and hence from (13^1) we obtain 

H 0 <1. (34) 

The models, in general, represent expanding, non-shearing and isotropic uni¬ 
verse. The models in the presence of bulk viscosity start expanding with a big 
bang at t = 0 when d = 0 and the expansion in the models decreases as time 
increases and the expansion stops at t = oo and t = — K 

4.2 Flat models 

For flat model Ricci scalar, R = 0 and hence in this case 

( k + at) 2 = [1 — a(^af 2 + kt + d)]. (35) 

The mass-density and pressure are given by 

8n P = ^2 + ai ) 2 “ X ] > ( 36 ) 

8tt( P - £6) = [(k + at) 2 - l] . (37) 

The expansion 6 is given 

3-i / [1 — a(iat 2 + kt + d)l 

9 = -V---• (38) 

Here we again consider two subcases: 
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4.2.1 Model I: solution for £ = £o 

When n = 0, Eq. (151)1) reduces to £ = £o = constant. Hence in this case Eq. 
(EZD, with the use of (1551) . leads to 



(39) 


4.2.2 Model II: solution for £ = £o p 

When n = l, Eq. (1551) reduces to £ = £o P- Hence in this case Eq. (EZl), with the 
use of EH) , reduces to 



(40) 


The models, in general, represent expanding, non-shearing and isotropic uni¬ 
verses. The models in the presence of bulk viscosity start expanding with a big 
bang at t = 0 when d = 0 and the expansion in the models decreases as time 
increases and the expansion stops at t = oo and a = 1. 

I Q 

5 Solution for b = — where a is constant 

In this case, on integrating, Eq. © gives the exact solution 


o 

where k and d are constants of integration. 



(41) 


Here we consider the two following cases: 


5.1 Non-flat models 


The mass-density, pressure and Ricci scalar are given by 



(42) 



(43) 



(44) 


The expression for kinematical parameter expansion 6 is given by 



( 45 ) 
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where S is given by m- 

For simplicity and realistic models of physical importance, we consider the 
following two subcases: 


5.1.1 Model I: solution for £ = £o 

When n = 0, Eq. (P2S1) reduces to £ = £o = constant. Hence in this case Eq. 
(l43l) . with the use of (l45l) . leads to 


247T& 1 2 1 

8nP = + 2^ ’ 


(,k + ^at * 2 ) 2 — 1 + 2 at(^at 3 + kt + d) 


■ (46) 


5.1.2 Model II: solution for £ = £o p 

When n = 1, Eq. (1^51) reduces to £ = £oP■ Hence in this case Eq. (H51) . with the 
use of TO . reduces to 


8np = 


{(fc + 4af 2 ) 2 - 1} 
S 3 


9^o (k + -at 2 )-S 


2a(^at 3 + kt + d) 
52 ' 


(47) 


5.1.3 Physical behaviour of the models 

From Eqs. (l42l) and (|44)l . it is observed that p > 0 and R > 0 for t > 

\J (2/a)(l — k), where k < 1. From ([45]) . we observe that 9 decreases as t 
increases. 


We find that shear cr = 0 in the models. Hence cr/9 = 0 which shows that 
the models are isotropic. At any intermediate time t = (2/a)(2 — k), R > 0 

implies that 

3 + ^(2 — fc)( 1 + k) + d^/2a{2 - k) > 0, (48) 

o 

where k <2. From Eq. ([48]) . it is evident that 


d\/2a(2 — k) > 0. 


(49) 


From Hubble’s parameter, H = ■|, we obtain a cubic equation in to 

3 3to 2 , 6 kt 0 6 

t 0 - ~rr- H-et ( fc _ dH °> = °- 

H 0 a aH o 


Solving Eq. © we obtain 


to = —! a = 2 kHo 2 and d = - 

H 0 3 Hi 


3 > 


0 


where to is an epoch time. Thus from Eq. © 

4£T 0 2 fc(2 - k) > 0. 


(50) 


(51) 


(52) 
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Hence Eq. USD, for k < 2, implies that 

H 0 > 0. (53) 

The models, in general, represent expanding, non-shearing and isotropic uni¬ 
verse. The models in the presence of bulk viscosity start expanding with a big 
bang at t = 0 when d = 0 and the expansion in the models decreases as time 
increases and the expansion stops at t = oo and t 2 = — 


5.2 Flat model 

For flat model Ricci scalar, R = 0 and hence in this case 
(,k + ^ai 2 ) 2 = [1 — at(-at 3 + kt + d)]. 
The mass-density and pressure are given by 


8 "’= s* 


8t t(p - £0) = 


(k + \at 2 ) 2 - 1 


S 2 


(k + ) ~ ^ 


The expansion 9 is given 


0 = 


3^/ [1 — at(^at 3 + kt + d)] 


S 


Here we again consider two subcases: 


(54) 


(55) 

(56) 


(57) 


5.2.1 Model I: solution for £ = £o 


When n = 0, Eq. (|26l) reduces to £ = £o = constant. Hence in 
([56]). with the use of ([57]) . leads to 


8irp = 


247t£ 0 


[1 - at(^at 3 + kt + d)\ + 


(k+^at 2 ) 2 - 


this case Eq. 
ll . (58) 


5.2.2 Model II: solution for £ = £o p 


When n = 1, Eq. (l26l) reduces to £ = £o p. Hence in this case Eq. (l56ll . with the 
use of ([57jl . reduces to 


8np = 


{(k + \at 2 ) 2 - 1} 

S 5 


9£o \/ [1 — at(-at 3 + kt + d)] + S 


(59) 


The models, in general, represent expanding, non-shearing and isotropic uni¬ 
verse. The energy conditions p > 0 and p > 0 are satisfied when sj (2/a)(l — k). 
The models in the presence of bulk viscosity start expanding with a big bang at 
t = 0 when d = 0 and the expansion in the models decreases as time increases 
and the expansion stops at t = oo and t 2 = — 
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6 Solution for b = — 4f-, where K is constant 

D 6 7 


In this case, on integrating, Eq. ([20]) gives the exact solution 

a , (a + 2I<t ) 3/2 
b ~ P + 3 K ’ 

where a and j3 are constants of integration. 


(60) 


Here we consider the two following cases: 


6.1 Non-flat models 


The mass-density, pressure and Ricci scalar are given by 


8lr P = 32 K a + 2Kt ) ~ !] > 


8tt( P - £ 0 ) = 


S 2 


1 - (a + 2Kt) - 


2 K 


R S 2 


{a + 2 Kt) - 1 + 


(a + 2Kt) 1 / 2 

I< 


{P- 


(a + 2 Kt) 2 
3 K 


a+ 2 Kt) 1 / 2 




The expansion 9 is calculated as 


= l {a + 2Kt)^. 


(61) 

(62) 

(63) 


(64) 


For simplicity and realistic models of physical importance, we consider the fol¬ 
lowing two subcases: 


6.1.1 Model I: solution for £ = £o 

When n = 0, Eq. (l26l) reduces to £ = £o = constant. Hence in this case Eq. 
®2D, with the use of (l 66 l) . leads to 


87rp 


24tt£ 0 

S 


(a+ 2 AT) 1 / 2 _J_ 


a + 2 Kt - 1 


2K r,o , (« + 2AT)l ' 
(a + 2AT) 1 / 2tP 3K 1 

(65)' 


6.1.2 Model II: solution for £ = £o P 

When n = 1, Eq. (J26]) reduces to £ = £o p. Hence in this case Eq. (l62l) . with the 
use of ((Ml) , reduces to 


8np 


(.a + 2Kt - 1) 
S* 


9(a + 2Kt) 1 / 2 - S 


2I< (a + 2Kt)i 

S 2 (a + 2Kt) 1 / 2XP 3K 

( 66 ) 
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6.1.3 Physical behaviour of the models 

From Eq. (ED, it is observed that p > 0 for t > (1 — a)/2I\, where a < 1. From 
ED, we observe that 6 decreases as t increases. 

We find that shear a = 0 in the models. Hence a/6 = 0 which shows that 
the models are isotropic in nature. 

From the Hubble’s parameter, H = )|, we obtain 

x ° 3 - 3K (jf 0 - p) = °> E) 

where xo = V a + 2Kto and to is an epoch time. For the solution of Eq. (ED 
let us assume xo = P and hence we have 

P 2 = 3 K (^ - l) . (68) 

From Eq. (l63l) . R > 0 implies that 

4x 3 - 3x + 3pK > 0, (69) 

where x = Vet + 2 Ixt. 

Solving Eq. (1M1) at any intermediate time, t = (3 — a)/2K, where a < 3, we 
obtain 

3 + IF > 0. (70) 

From xo = V a + 2/Fto = P, we have 

a + 2A ^=^ 

which reduces to 

P 2 = 3. (71) 

From Eqs. (l68l) and (I7T1) . we obtain 

K =^—. (72) 

Ho 1 

From Eqs. GOD and ED, we have 

H 0 < 1.5. (73) 

From ED in order that I\ > 0, Hq < 1. Hence the value of Hubble’s constant 

is less than unity. 
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6.2 Flat model 


For flat model, Ricci scalar R = 0 gives at any intermediate time, t = (3 — 
a.\)/2k\, where ay < 3, 


3 H- k\ — 0. 

(74) 

The above equation with the use of Eq. (l72l) gives 


H 0 = 1.5 

(75) 

The mass-density and pressure are given by 


Snp = [(ai + 2 kit) - 1 ] , 

(76) 

8 t r(p - £6) = [(ay + 2 kit) - 1 ] , 

(77) 

where 

O O , Oi + 2 M ) 3/2 

S = Pl+ 3 Ay 

(78) 

and k\ , ol\ and /3i are constants. In non-flat model, Eq. 
and from (1751) using H 0 = 1.5, we have from (|72|). k\ = - 
positive for 

t> [a i - (3/?i fc 2 ) 2/3 ) /2fci, 

(1711) gives pi = V 3 
-3. Accordingly S is 

where &2 = 3 such that (3/?ifc2) 2 / 3 < ai <3. 



7 Conclusions 

In this paper we have described a new class of LRS Bianchi type I cosmological 
models with a bulk viscous fluid as the source of matter by applying variable 
deceleration parameter. Generally, the models are expanding, non-shearing and 
isotropic in nature. The models in the presence of bulk viscosity start expanding 
with a big bang at t = 0 when d = 0 and the expansion in the models decreases 
as time increases and the expansion stops at t = oo and t = —K The study of 
the results of the three deceleration parameter models of the universe showed 
that the Hubble’s constant is less than 1.00 for the non-flat model, and it is 
1.5 for the flat model. These mathematical results are of the physical interest. 
The scale factor is not obtained to be linearly related to time as the case with 
supernova cosmology model [40]. If £ = 0 is set in the solutions obtained in 
this paper, we get the solutions obtained by Paul [39]. But in his paper [39]. 
there are errors in the the beginning equations (1151) and (ll(il) which propagate 
through out the paper and affect all the results. 

The coefficient of bulk viscosity is taken to be a power function of mass 
density. The effect of the bulk viscosity is to produce a change in perfect fluid 
and hence exhibit essential influence on the character of the solution. Murphy 
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[26] has studied perfect fluid cosmological models with bulk viscosity and ob¬ 
tained that the big bang singularity may be avoided in the finite past. We also 
observe here that Murphy’s condition [7] about the absence of a big bang type 
singularity in the finite past in models with bulk viscous fluid, in general, not 
true. 
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